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ABSTRACT 

We compare the H I column densities from Lya absorption to the metal column densities from X-ray ab- 
sorption in 7-ray burst (GRB) afterglows. Eight bursts of seventeen entering the sample, have significant 
extragalactic X-ray absorptions. Of these we find a range in metallicities from sub-solar to a few hundred times 
the solar value. There is a complete lack of correlation between the X-ray and optical column densities. This, 
and the large range and extreme values of these 'metallicities', suggest that the column densities derived by 
one or both methods are not a reliable indication of the true total column densities towards GRBs. Ionisation 
of the GRB's gas cloud to large distances along the line of sight seems the most likely cause. From the lower 
limit on the total column density and the UV luminosity of the GRBs, we derive a maximum distance to the 
majority of the gas surrounding GRBs of ~ 3 pc, suggesting that the gas probed by optical afterglow spectra is 
not the cloud in which the burst occurs. This is an encouraging result for the use of GRB optical afterglows as 
probes of the interstellar medium (ISM) in their host galaxies, as the ISM observed is less likely to be strongly 
affected by the GRB or its progenitor. 

Subject headings: gamma rays: bursts - X-rays: general - X-rays: galaxies - quasars: absorption lines 



1. INTRODUCTION 

It is known from studies of neutral hydrogen absorption in 
7-ray burst (GRB) afterglows (a technique pioneered half a 
decade ago, Jensen et al. 2001; Fynbo et al. 2001; Hjorth 
et al. 2003), that they probe much higher column densities 
than QSO sightlines (Vreeswijk et al. 2004), almost certainly 
because GRBs probe the star-forming environments in which 
they are born (Jakobsson et al. 2006a). The highest column 
density GRB damped Lya absorbers (DLAs) are higher than 
any known QSO-DLA (Jakobsson et al. 2006a; Watson et al. 
2006a; Starling et al. 2005; Chen et al. 2005), in spite of the 
much larger QSO-DLA sample. 

Measurements of the H I column density (N(H I)) are 
needed to determine abundances in GRB environments where 
the metal column densities are known either from absorption 
lines of non-refractory elements such as Zn, or from the soft 
X-ray absorption which is dominated by a-chain elements in 
gas and solid phases. The optical observations are compli- 
cated by the need to assume that certain elements are not de- 
pleted by dust, and the need for high-resolution spectra. The 
X-ray measurements do not suffer from these problems, but 
do not provide a strong redshift constraint on the absorber 
(e.g. Watson et al. 2002; Stratta et al. 2004; Gendre et al. 
2006). However using the X-ray metal column densities with 
Lya columns has been extremely difficult because of the con- 
flicting redshift requirements and has not yet been done in any 
object at z > 0.5 other than a GRB. To detect hydrogen Lya 
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1216 A with ground-based instruments we need z > 2 to shift 
the line into the sensitivity range of optical/UV spectrographs. 
But the soft X-ray absorption moves further and further out 
of the bandpass as we go to higher redshift, such that once 
the neutral oxygen edge at 0.52 keV is no longer detectable, 
there is typically a degeneracy between the redshift and the 
total absorption column that scales roughly as A/h oc (1 +z) 2A . 
Thus we need a high X-ray column density (greater than a few 
x 10 21 ) for it to be detectable at z > 2. 

With the advent of the Swift satellite and its highly success- 
ful X-ray telescope (XRT), very early, largely complete, and 
high flux observations of GRB X-ray afterglows have been 
obtained for about ninety GRBs every year. This means that 
the sample of GRB afterglow X-ray absorption column den- 
sities has increased enormously, such that we now have limits 
on the extragalactic X-ray absorption for almost every new 
GRB discovered and we finally have a few GRB damped Lya 
systems (GRB-DLAs) where we can compare with signif- 
icant detections of extragalactic soft X-ray absorption: see 
GRB 050401 (Watson et al. 2006a), GRB 050505 (Hurkett 
et al. 2006), and GRB 050730 (Starling et al. 2005). 

In this Letter we present the sample of GRB-DLAs with 
significant extragalactic X-ray absorptions and compare the 
column densities obtained with each method. Uncertain- 
ties quoted are at the 68% confidence level for one interest- 
ing parameter unless otherwise stated. A cosmology where 
Ho = 70km s" 1 Mpc" 1 , fl\ = 0.7 and £l m = 0.3 is assumed 
throughout. 

2. SAMPLE AND DATA REDUCTION 

Swift GRBs at z > 2 were examined by Jakobsson et al. 
(2006a) to determine the properties of GRB-DLAs. From 
that sample bursts with well-constrained H I columns were 
selected. We analysed XRT data from these bursts where they 
were believed to have X-ray absorption in excess of Galactic 
from the literature or the Swift data table. 6 Analysis was done 
in a standard way. Spectra were fit with a power-law with 

6 http://swift.gsfc.nasa.gov/docs/swift/archive/grb_table/ 
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absorption fixed at the Galactic value and a second absorber 
at the host redshift. Where the count rate was high, regions 
with the PSF-core excised were used for spectral analysis to 
account for pile-up. 

To ensure that only statistically significant excess absorp- 
tion columns were selected, two conservative criteria were 
applied. 1) GRBs at low Galactic latitude (\b\ < 20°) were 
excluded. 2) The excess column density had to be detected 
above the 99.7% confidence level when the z = absorber 
was fixed at 1 x 10 20 cm" 2 or 20% above the Galactic value, 
whichever was larger (see Elvis et al. 1986, 1989; Dickey & 
Lockman 1990). Galactic column densities obtained from the 
newer H I survey of Kalberla et al. (2005) were consistent 
with those obtained from Dickey & Lockman (1990) within 
~ 3 x 10 19 cm -2 . At early times an absorbed power-law model 
fit occasionally yields excess absorption not present at later 
times (e.g. Starling et al. 2005; Gendre et al. 2007; Campana 
et al. 2007). Butler & Kocevski (2006) have examined the 
reasons behind this and we discuss it in § 4. Because of this 
effect we adopt the lower value of absorption where inconsis- 
tent values were found at early and later times. GRBs 050730, 
050820A, and 060714 are the bursts affected by this, and we 
have adopted absorption values from the later PC mode spec- 
tra for these bursts. 

The spectra of GRBs 050401, 050505, 050730 and 050904 
have been presented elsewhere (Watson et al. 2006a; De 
Pasquale et al. 2006; Hurkett et al. 2006; Starling et al. 2005; 
Watson et al. 2006b; Gendre et al. 2007; Campana et al. 2007). 
For the bursts with persistent excess absorptions published 
previously— 050319 (Campana et al. 2006), 050401, 050505, 
and 050904 — we obtain somewhat different values. This may 
be due to the more recent calibration used in this analysis. 
However, in each case, the authors claim a detection of ex- 
cess absorption and adopting those values do not change our 
conclusions. 

3. RESULTS 

Seventeen GRBs enter our sample as set out in Table 1. 
The comparison of X-ray and and hydrogen Lya absorption 
is shown in Fig. 1. If the metallicities of GRB formation sites 
are low (Le Floc'h et al. 2003; Fynbo et al. 2003; Fruchter 
et al. 2006; Modjaz et al. 2007; Kewley et al. 2007), we would 
expect to make no significant detections of excess absorption 
in this sample. Instead, eight GRBs have significant excess 
column densities: GRBs 050319, 050401, 050505, 050904, 
060210, 060607A, 060714, and 060926 (see Table 1 and 
Fig. 1). 

The best limits that can be placed on the column density 
from Swift-XRT data above z = 2 is roughly logA_i < 21.5. 
Fig. 1 is therefore not populated below 21.5 in X-ray col- 
umn. Such a limit could only constrain the metal abundance 
to be less than solar in the highest H I column density cases. 
Upper limits obtained in the optical/UV are more interest- 
ing: GRB 060607A has a tightly constrained, low H I col- 
umn density, but large X-ray absorption. (The same is true for 
GRB 060124, but it does not enter our sample at b = 17°.) 

There is evidence of a spread in the metallicities of GRB 
hosts, though with most fairly close to ~ 10% of the solar 
value (see Fynbo et al. 2006). If the metallicities of GRB en- 
vironments were all approximately similar we would expect 
some correlation between the optical and X-ray column den- 
sity determinations. It is clear from Fig. 1 that there is no such 
correlation and that the optical-to-X-ray ratios span a range 
of at least a few hundred. Furthermore, the metal abundances 



TABLE 1 

UV AND X-RAY ESTIMATES OF THE ABSORBING COLUMN DENSITIES IN 

GRBs 



GRB 


z a 


logjV H i fl 


Total /V H 

at 7, = 


Gal. N H 

(10 20 cm- 2 ) 


A'h atz 


050319 


3.24 


20.9 




1.1 


76+f> 

— 18 


050401 


2.90 


22.6 


16.9 ±0.8 


4.8 


196±13 


050505 


4.27 


22.1 


9.8 ±1.0 


2.0 


130 ±36 


050730 


3.97 


22.1 


3 


3.05 


< 100 


050820A 


2.61 


21.1 


4.0 ±1.0 


4.71 


< 178 


050904 


6.30 


21.3 


9.1 ±0.6 


4.9 


380 ±50 


050908 


3.34 


19.2 


2.2±1.5 


2.1 


< 122 


050922C 


2.20 


21.6 


6.5* 


5.75 


<55 r 


060206 


4.05 


20.9 


5 2 +lA 


0.9 


<280 


060210 


3.91 


21.7 


15.2±0.5 


8.5 


169 ±18 


060522 


5.11 


20.5 


4.6* 


4.83 


170 c 


060526 


3.22 


20.0 


6* 


5.51 


<84 r 


060607A 


3.08 


< 19.5 


5.6 ±0.4 


2.7 


55±7 


060707 


3.43 


21.0 


1.8* 


1.76 


<90 r 


060714 


2.71 


21.8 


1 3 6 +2 -° 


6.7 


94+27 


060906 


3.69 


21.9 


9.66* 


9.66 


< 100 c 


060926 


3.21 


22.7 


20.5^ 


7.31 


250+110 



" See Jakobsson et al. (2006a) and references therein. 
* From http://swift.gsfc.nasa.gov/docs/swift/archive/grb_table/ 
c Estimated using a redshift-corrected, 3 X 10 20 uncertainty in the X-ray col- 
umn and a 1 X 10 20 uncertainty in the Galactic column density. 

derived by this comparison indicate metallicities far above the 
solar values. The extreme abundances and the lack of correla- 
tion between the optical and X-ray values are disturbing and 
strongly imply that the X-ray absorption and the Lya absorp- 
tion are sampling different environments. These data immedi- 
ately show that for GRBs the X-ray-to-optical absorption ratio 
is not a useful measure of the metallicity. 

3.1. Reliability of the measurements 

Metal abundances have been determined for intervening ex- 
tragalactic absorbers in a few radio-loud AGN using X-ray 
and Lya absorption. Those values are significantly below so- 
lar (see Turnsheket al. 2003). The metallicities are consistent 
with extinction measures and expectations for these objects. 
For the Crab pulsar, this technique yields values consistent 
with solar metallicity (Sollerman et al. 2000; Willingale et al. 
2001). Both these results are plotted for comparison on Fig. 1 
and indicate that in general the comparison of column densi- 
ties from X-ray and Lya is valid for both Galactic and extra- 
galactic sources. 

GRB afterglows exhibit curvature close to the cooling fre- 
quency, which may be at X-ray wavelengths at certain times 
and may mimic absorption. Butler & Kocevski (2006) have 
shown that in many cases where the soft X-ray absorption 
apparently changes rapidly in the first few hundred second 
after the burst, this may be better explained by curvature of 
the continuum than ionisation of the absorber. This seems 
a reasonable explanation in the early phases of some bursts. 
Such early curvature however, should not affect the results 
presented here, as we have used absorption values from later 
spectra where early (high) values of the absorption are incon- 
sistent with later data. Such an approach is conservative, in 
the sense that it leads to lower column densities. In general, 
spectral curvature is not responsible for the detection of ex- 
cess absorption. Most of the bursts in question have early WT 
and PC mode data that agree on the absorption to within la 
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FIG. 1 . — Neutral hydrogen column density (from the UV) as a function of 
X-ray equivalent hydrogen column density. The ratio of these values provides 
an estimate of the metal abundance of the absorbing medium. The clear lack 
of correlation between these column densities seems initially to indicate an 
extremely strong variation in the metal enrichment of the absorbing clouds 
near GRBs. These values are so extreme however, and many so much greater 
than solar, that it seems likely that there are strong systematic uncertainties 
related to either or both of the absorbing column density estimates from X- 
rays and the UV, possibly because of ionisation of most of the line-of-sight 
hydrogen. Data from absorption toward a Galactic source (the Crab pulsar) 
and a high column density DLA (foreground to the blazar AO 0235+164) are 
shown for comparison. Where the uncertainties on the H I column is not in 
the literature a value of 0. 1 has been plotted. 

despite large spectral changes. 

Given that the technique is reliable, could it be that the 
Galactic column densities are underestimated? This explana- 
tion can be ruled out. First of all we have excluded any excess 
column density that is not significant above the 90% varia- 
tions of the Galactic columns (Elvis et al. 1986, 1989; Dickey 
& Lockman 1990). Furthermore solar abundances were used 
to convert the Galactic 21cm column densities (Dickey & 
Lockman 1990), resulting in a high (conservative) estimate 
of the foreground Galactic column, higher than is often ar- 
gued for the Galactic gas by e.g. Wilms et al. (2000). A large 
underestimate of the Galactic column density can also be ex- 
cluded on the basis of observations of blazars with BeppoSAX 
(Donato et al. 2005) that show no substantial absorption com- 
pared to the values of Dickey & Lockman (1990). Of nearly 
90 blazars examined with BeppoSAX only 17 showed ex- 
cess absorption in the observed 0.1-50keV spectra and even 
these excesses are small, ~ 1 x 10 20 cm" 2 . Finally, if the ex- 



cess columns were due to Galactic foreground absorption, we 
would expect a systematic rise in observed column densities 
with redshift, an effect we do not observe when we combine 
with the lower-redshift values of Campana et al. (2006). 

4. DISCUSSION 

It is only with the success of Swift that we have been able to 
obtain Lya and X-ray absorption measures in the same bursts. 
Given the expectations for neutral hydrogen column densi- 
ties surrounding GRBs (Reichart & Price 2002), few signif- 
icant detections of excess absorption in the X-ray afterglows 
of GRBs should have been expected above z > 2. 

4.1. Intervening absorbers 

It is expected that high redshift sources (and GRBs have 
a very high mean redshift, Jakobsson et al. 2006b) will of- 
ten have intervening line of sight absorbers. The X-ray ab- 
sorption does not tell us about the redshift of the absorber 
and some absorption could in principle be related to low-z 
systems. The observed absorption in X-rays drops substan- 
tially as the absorber is moved to higher redshift, so for this 
to be an effective explanation any absorber would have to be 
at relatively low redshift. However, very large column densi- 
ties of metals would routinely be required along most lines 
of sight at low redshifts. Recently, Prochter et al. (2006) 
have shown that GRBs have intervening large column density 
(Wa > 1 A) Mg II absorbers about four times more frequently 
than QSOs. And it is a suggestive coincidence that GRBs also 
have significant excess soft X-ray absorptions. However, to 
explain the results observed here with intervening absorbers 
is not possible. This would require an absorber with typ- 
ically log(A^n) — 22.1 at z = 0.3 (assuming Z/Z = 0.1) in 
half of the GRB sightlines — an extremely large value. As- 
suming that an Mg II absorber with Wa ~ 1 - 3 A corresponds 
to log(/V H ) < 21 at this metallicity (Meiring et al. 2007), the 
required absorption is roughly two orders of magnitude above 
that provided by the Mg II absorbers seen in GRBs at z < 0.5 
(Prochter et al. 2006). To reinforce this point, a similar ex- 
cess of large equivalent width Mg II absorbers is also found 
in the foregrounds of blazars (Stocke & Rector 1997), which, 
as noted above, do not have large X-ray absorptions (Donato 
et al. 2005). This is an important difference in that it seems 
likely that GRBs and blazars sample similar column densities 
of intervening absorbers, but GRBs typically have large soft 
X-ray absorptions and blazars do not, showing that the large 
absorptions are indeed intrinsic to the GRB hosts. 

It has been suggested (Jakobsson et al. 2006a) that the lack 
of very high column density and apparent overabundance of 
low column density sources in comparison to model predic- 
tions may be due to ionisation of the hydrogen near the GRB. 
This hypothesis explains the results found here; the hydrogen 
may be ionised to large distances by the early afterglow, but 
the metals, while also ionised, will not be stripped to such an 
extent that they cease substantially to absorb the soft X-rays. 
This hypothesis is also consistent with the lack of expected 
Wolf-Rayet features in most GRB optical afterglows (Chen 
et al. 2006). 

4.2. Size of the gas cloud 

It is possible to limit the size of the gas cloud where most 
of the hydrogen has been ionised by the GRB. A maximum 
distance to the gas is found by assuming the X-ray column 
density represents the total column density (it cannot obvi- 
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ously be smaller than this, and if it is larger, then the ionisa- 
tion radius will be smaller) and that each ionising UV photon 
is intercepted by a hydrogen atom (if this is not the case, the 
radius will again, necessarily be smaller). Since we know the 
total UV fluence for most of these bursts is a small fraction 
of the prompt energy, we can derive a maximum distance to 
which the GRB could have ionised the gas. In all cases, this 
is less than ~ 3 pc. Interestingly, if the gas cloud is already 
ionised by the massive stars in the region, the typical densi- 
ties of H II regions (< 10 4 cm" 3 ) yield sizes which are also at 
most a few parsecs in radius. Smaller radii require higher den- 
sities. Radii an order of magnitude greater than this make the 
total masses unfeasibly large. This radius is similar to that ob- 
tained by Campana et al. (2007) based on the assumption that 
the apparent decreasing absorption detected in GRB 050904 
is due to ionisation by the GRB. 

Observations of UV absorption lines have allowed dis- 
tances to low-ionisation gas to be derived that are far larger 
than the maximum distance to the majority of the gas col- 
umn we derive above. Prochaska et al. (2006) find > 50 pc 
for Mg I in GRB 0501 11, and Vreeswijk et al. (2006) find 
1.7±0.2kpc for Ni II and Fe II in GRB 060418. Gas with 
high ionisation states closer to the GRB has been detected in 
the optical (D'Elia et al. 2006). Though 1.7 kpc seems large 
to be the typical distance scale between clouds in most bursts, 
taken together, those results support our conclusion in this pa- 
per, that the optical/UV spectra are probing gas that is outside 
the pc-scale structure in which the GRB explodes. 

5. CONCLUSIONS 

We have examined the X-ray afterglow spectra of the sam- 
ple of Swift GRBs at z > 2 as defined by Jakobsson et al. 
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(2006a) and compared their extragalactic soft X-ray absorb- 
ing column densities to the host H I column densities. We 
find no correlation between the column density values, and 
the X-ray absorptions often far exceed the H I column densi- 
ties. The most likely explanation for this discrepancy is that 
GRBs substantially ionise the hydrogen in their immediate 
environments (at least along the line of sight). This may be 
consistent with the sublimation of the dust that is expected (by 
analogy with the local universe) to be associated with large X- 
ray column densities, though it does not provide evidence in 
favour of dust sublimation. It seems likely that significant in- 
formation about the immediate circumburst environment may 
be gleaned from high-resolution X-ray spectroscopy of the 
very early afterglow, probably revealing the complex absorp- 
tion pattern associated with metal-rich warm absorbers, as ob- 
served in some AGN. Finally, we find the maximum distance 
for the bulk of the absorbing material surrounding the GRBs 
to be ~ 3 pc for these bursts, similar to the sizes of compact 
H II regions. Since the distances inferred from optical obser- 
vations are typically larger than this, it implies that the opti- 
cal/UV spectroscopy of GRB afterglows typically probes en- 
vironments that are little affected by the progenitor, the explo- 
sion or the bright afterglow. This is particularly encouraging 
for the use of GRBs as probes of star-forming regions at high 
redshift. 
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